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Section  I.  INTRODUCTION 


The  mallaunch  of  free  rockets  continues  as  a significant  problem. 

Each  new  system,  i.e.,  each  new  launcher /rocket  combination,  possesses 
unique  characteristics  which  must  be  analyzed  and  understood  before  the 
performance  of  the  system  can  be  predicted  or  optimized.  The  analyses 
which  are  conducted  are  usually  both  analytical  and  experimental.  Simple 
physical  models  of  the  actual  system  are  usually  devised  and  mathematical 
theories  of  how  the  system  will  behave  are  constructed  on  the  bases  of 
such  physical  models.  The  behavior  of  a launcher /rocket  system  predicted 
by  such  a theory  is  normally  compared  with  experimental  data.  By  doing 
this,  improvements  in  both  the  theory  and  the  design  of  the  system  can 
usually  be  made.  Furthermore,  once  a theory  which  predicts  the  behavior 
of  a system  or  class  of  systems  is  devised,  analytical  experiments  using 
the  theory  often  give  impetus  to  the  development  of  radically  new  systems. 

The  work  reported  herein  represents  a first  step  in  developing  a 
theory  for  the  dynamical  behavior  of  free  rocket  launch  systems.  As 
such,  the  results  of  this  study  should  be  useful  in  the  ways  previously 
mentioned.  In  particular,  this  study  appears  to  be  the  first  in  which 
the  motion  of  rocket/launcher  systems  during  the  detent,  guidance,  tip- 
off,  and  free- flight  phases  of  the  complete  flight  is  considered. 

Studies  of  various  aspects  of  free  rocket  launch  dynamics  have 
already  been  made  [1,2,3].  Because  of  the  advances  in  computational 
methods  and  simulation  of  dynamic  systems  in  general  it  is  important 
that  such  a study  be  conducted.  This  report  is  a contribution  in  this 
regard . 

In  Section  II,  a simple  physical  model  for  a single-round 
launcher/rocket  system  is  described.  This  model  includes  the  effects  of 
many  factors  which  contribute  to  mallaunch  of  free  rockets  and  possesses 
enough  flexibility  to  allow  for  addition  of  the  effects  of  other  factors. 

Differential  equations,  the  solution  of  which  provides  a description 
of  the  motion  of  the  system,  are  derived  in  Section  III  using  the 
Eulerian  vector  approach.  The  vector  equations  which  are  obtained  are 
converted  to  matrix  equations  amenable  to  digital  solution  by  using 
modern  matrix/vector  algebra.  As  a result  of  such  conversion,  the  burden 
of  performing  much  of  the  algebra  involved  i*n  computing  derivatives  of 
the  dependent  variables  may  be  placed  on  the  digital  computer. 

Section  IV  provides  descriptions  of  some  of  the  methods  used  in 
implementing  the  digital  program.  In  particular,  the  manner  in  which 
spin  torque  and  thrust  profiles  are  generated,  matrix  algebra  is  per- 
formed and  the  integration  of  the  equations  is  carried  out  are  explained. 

Typical  results  obtained  using  the  digital  computer  program  (Fortran 
IV  computer  code)  are  presented  and  discussed  in  Section  V and  a summary 
and  conclusions  based  on  the  study  are  presented  in  Section  VI.  The 
appendices  are  descriptions  of  the  input  needed  for  the  program  and  a 
listing  of  the  computer  code. 
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Section  II.  PHYSICAL  MODEL  FOR  THE  SYSTEM 


a.  General  Comments 


The  physical  model  adopted  for  the  current  study  is  described 
in  this  section.  The  model  is  fairly  general  in  that  it  allows  for 
rotational  motions  of  an  arbitrary,  rigid-body  launcher  and  arbitrary 
(six  degrees  of  freedom  when  not  constrained  by  the  launcher  and  five 
degrees  of  freedom  when  constrained)  motion  of  the  rocket  which  is 
modeled  as  a constant  mass  rigid  body.  The  model  includes  the  effects 
of: 

1)  Thrust  malalignment. 

2)  Large  initial  launch  angle  (elevation) . 

3)  Mass  unbalance  (both  static  and  dynamic)  of  the  rocket. 

4)  Flexibility  of  the  launcher /rocket  interface. 

5)  The  detent  force. 

6)  Gravity  tipoff. 

7)  Variable  spin  and  thrust  programs. 

8)  Friction. 

Figure  1 is  a graphical  depiction  of  the  complete  model.  In  that 
figure  is  shown  a rocket  of  mass  m,  which  rests  upon  a launcher  of  mass 
M.  The  rocket  is  connected  to  the  launcher  by  mechanical  devices  which 
rigidly  constrain  the  rocket1 s relative  motion  along  the  prescribed  line 
of  flight  (through  a detent  device  which  is  assumed  to  act  at  the  rear 
rocket  support)  and  flexibly  constrain  relative  pitch  and  yaw  motions  of 
the  rocket.  The  system  of  two  rigid  bodies  adopted  for  the  model  has  up 
to  nine  degrees  of  freedom,  with  eight  during  the  detent  phase. 

In  the  following  sections,  the  launcher  and  rocket  models  are 
treated  separately  and  relevant  definitions  are  given. 

b.  Launcher  Model 


The  launcher  model  may  be  represented  by  the  single  rigid  body 

shown  in  Figure  2.  This  body  has  three  degrees  of  rotational  freedom 

about  the  fixed  point  0 which,  in  general,  does  not  coincide  with  the 

launcher's  center  of  mass  C . The  coordinate  system*  OXYZ , which  is 

L 

inertial  and  has  its  Z-axis  directed  vertically  down  and  its  X-axis 


*All  coordinate  systems  used  or  referred  to  herein  are  dextral, 
orthogonal  systems. 
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Figure  1.  System  model. 
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Figure  2.  Launcher  physical  model. 


directed  parallel  to  the  horizontal  projection  of  the  desired  flight 

path  of  the  rocket,  and  the  C x y z system,  which  is  fixed  in  the 

L L L L 

launcher,  are  also  shown  in  Figure  2.  The  stiffness  and  damping  of  the 
launcher  are  assumed  to  be  concentrated  at  the  point  0,  exerting  a 
torque,  , about  the  point.  The  mass  distribution  of  the  launcher  may 

be  arbitrary,  but  is  assumed  to  be  known. 
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The  springs,  which  are  shown  connecting  a rocket  to  the  launcher, 

are  intended  to  model  the  flexibility  characteristics  of  the  parts  of  the 

launcher  and  the  rocket  which  couple  the  two,  i.e.,  the  launcher /rocket 

interface.  After  detent  release  the  springs  translate  parallel  to  the 

x-axis  with  the  rocket  and  exert  forces  in  the  y - and  z -directions  at 
Li  Li  Li 

the  points  f and  a.  Also,  frictional  forces  resulting  from  sliding  con- 
tact of  the  shoes  (or  other  connecting  devices)  are  assumed  to  act, 
opposing  the  sliding  motion.  Finally,  gravitational  effects  predicated 
on  a flat  earth  model  are  included. 

In  Figure  3,  the  vector  R locates  the  launchers  center  of  mass  and 

R and  S together  locate  the  point  a,  which  lies  on  the  centerline 

— — a 

(geometric)  of  the  missile  and  the  y z plane  which  contains  the  aft 

Li  Li 

shoes.  Furthermore,  Mg  is  the  launcher  weight  vector, 
c . Rocket  Model 


The  rocket  model  which  was  adopted  for  this  study  is  shown  in 
Figure  3.  The  rocket  is  assumed  to  be  a single  rigid  body  of  constant 
mass,  m,  and  the  mass  distribution  of  the  rocket  is  assumed  to  be 
arbitrary.  The  forces  which  are  considered  to  be  acting  on  the  missile 
are  the  thrust,  the  rocket fs  weight,  forces  transmitted  by  the  springs 
which  connect  the  rocket  to  the  launcher,  frictional  forces,  the  detent 
force  (which  acts  at  a),  and  external  forces,  such  as  pressure  forces  due 
to  rocket  plume  impingement  on  the  launcher  and  aerodynamic  forces.  The 
torques  which  act  on  the  rocket  include  those  due  to  the  above  forces  and 
the  spin  torque  generated  by  a spin  motor  or  other  device  and  frictional 
torques  in  the  spin  bearings  (or  tube  reactions  in  a tube-type  launcher).* 

In  Figure  3,  the  vector  u + S + R locates  the  rocket's  center  of 

**  “3  “ 

mass,  which  does  not  necessarily  lie  on  the  rocket's  geometric  axis  of 

symmetry.  The  C x y z coordinate  system  is  fixed  in  the  rocket,  and 
m mm  m 

the  vectors  £ ^ and  £ ^ will  be  used  in  following  sections  in  the  deriva- 
tion of  equations  of  motion. 


*A  nonspinning  rocket  model  with  the  shoes  modeled  as  springs  has  been 
briefly  studied  in  Reference  4. 
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Figure  3.  Rocket  mathematical  model. 


Section  III.  EQUATIONS  OF  MOTION 


a.  General  Comments 


All  derivations  of  equations  of  motion  for  dynamical  systems 
in  which  relativistic  effects  may  be  ignored  are  based  on  Newton's  laws 
of  motion.  Sometimes  these  laws  are  applied  directly,  an  approach  which 
we  term  "Eulerian."  Other  times  the  dynamicist  may  use  the  formalism  of 
the  Lagrangian  and  Hamiltonian*  approaches  in  applying  Newton's  laws 
somewhat  indirectly.  The  choice  of  approach,  or  method  of  attack,  is 
somewhat  arbitrary  depending  on  the  type  of  constraints,  the  purpose 
for  which  the  equations  are  being  obtained,  the  prejudices  of  the  analyst, 
etc . 


When  the  constraints  on  a fairly  complex  dynamical  system  are  holo- 
nomic,  Lagrange's  equations  are  often  used  because  this  allows  the  con- 
straints to  be  more  or  less  automatically  satisfied  and  because,  once  the 
kinetic  energy  and  generalized  forces  have  been  formulated,  the  rest  of 
the  derivation  is  essentially  mechanical. 

The  Lagrangian  approach  has  drawbacks,  however.  For  example,  the 
second-order  differential  equations  which  are  obtained  thereby  are  often 
intricately  coupled  and  very  complicated  and  the  algebra  involved  in  the 
derivation,  although  mostly  mechanical,  may  be  extensive.  Furthermore, 
if  the  system  model  is  changed  in  some  manner,  it  may  be  necessary  to 
repeat  a large  portion  of  the  derivation. 

The  Eulerian  approach's  major  drawback  is  probably  the  difficulty 
of  including  the  effects  of  constraints  into  the  equations.  However, 
the  Eulerian  approach  is  straightforward  and,  in  the  author's  opinion, 
offers  more  flexibility  than  the  Lagrangian  approach,  especially  when 
modern  vector /matrix  methods  are  used  to  convert  vector  equations  of 
motion  to  matrix  equations  which  are  to  be  solved  numerically  on  a digi- 
tal computer.  Hence,  the  Eulerian  approach  was  adopted  for  the  deriva- 
tion which  follows  in  Subsections  III.c  through  Ill.f. 

As  a prelude  to  the  derivation,  some  comments  on  vector /matrix  nota- 
tion are  made  in  Subsection  Ill.b.  Next,  in  Subsection  III.c,  Newton's 
laws  are  applied  to  obtain  the  governing  dynamical  equations  in  first- 
order  form.  The  kinematical  aspects  of  the  launcher /rocket  system  are 
treated  in  Subsection  Ill.d,  while  the  effects  of  constraining  the 
rocket's  motion  during  the  detent  phase  of  launch  are  discussed  in  Sub- 
section Ill.e  and  explicit  expressions  for  the  external  forces  and  torques 
acting  on  the  rocket  and  the  system  are  determined  in  Subsection  Ill.f. 


^Usually  the  Hamiltonian  approach  is  reserved  for  use  when  the  system 
is  conservative. 
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b . Vector /Matrix  Notation 


In  deriving  equations  of  motion  for  complex  dynamical  systems, 
vectors,  although  not  completely  indispensable,  are  very  nearly  so. 

Vector  differential  equations  cannot,  of  course,  be  directly  programmed 
for  numerical  solution;  however,  one  may  transform  such  vector  equations 
in  equivalent  matrix  equations  in  which  matrix  operations  isomorphic  to 
vector  operations,  such  as  multiplying  a row  matrix  times  a column  matrix 
to  get  the  result  of  a vector  dot  product,  are  used.  Moreover,  by  using 
matrix  operators  isomorphic  to  vector  operators,  much  algebra  may  be 
avoided  by  the  analyst  and  passed  on  to  the  computer. 

As  an  example,  we  consider  the  vector  equation  of  motion  for  a point 
mass  m moving  under  the  influence  of  an  external  force  F in  a coordinate 
system  Oerj£  which  has  a fixed  origin  0,  but  is  rotating  with  angular 
velocity  w and  angular  acceleration  u)  as  shown  in  Figure  4. 


m 


Figure  4.  Geometry  for  the  point  mass  motion  example. 


The  equation  of  motion  may  be  written  in  the  vector  form, 

mr  = m[°v  + u)xr  + 2u)xr+D  x(to  x r)  ] = F , (1) 
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where 


oo 

r = 

i eL  + ijS2  + ?e3 

(2) 

OJ  = 

Wlel  + ^2  + W3e3  * 

(3) 

0 

r = 

• • 

t,  A • A t A 

5 e1  + qe2  + ^e3 

(4) 

rotate  with  the  0£t)£  system. 

Equation  (1)  holds 

regardless  of  the  system  in  which  F is  expressed.  For  example,  we  may 
write  F as 


F * Fx1  + V + FZK 


(5) 


A A A 

where  I,  J,  and  K are  nonrotating  unit  vectors  associated  with  the  iner- 
tial system  OXYZ  (Figure  4) . 


To  write  Equation  (1)  in  an  equivalent  matrix  form,  we  must  first 
express  each  terra  in  the  same  basis;  i.e.,  each  vector  term  must  be 
expressed  in  terms  of  the  same  unit  vectors.  If  C is  the  direction 
cosine  matrix  such  that 


then 


F = F1  e]_  + F2  e2  + F3  e3 


where 


(6) 


(7) 


(8) 


11 


Using  the  same  notation  for  column  matrices  and  the  vector  counter- 
parts, e.g.. 


we  may  write 


oo  'o 

m(r  + u r + 


2w  r + Qjur) 


C F 


(9) 


(10) 


where  the  tilde  above  a vector  denotes  a certain  skew- symmetric  matrix 
formed  from  its  components.  To  illustrate  the  tilde  notation,  we 
consider  [5] 


U) 


0 -w3  u>2 

0 - 

- w2  w 0 


(ID 


Upon  multiplying  u)  times  the  column  matrix  r,  we  find  that  the  elements 
of  the  product  column  vector  are  the  components  of  wxr  in  the  (£,  j,  £) 
basis . 

T o 

Returning  to  Equation  (10),  solving  it  for  jr  and  defining  (u  v w)  =_r 
(where  T denotes  the  transpose) , we  obtain 


(u  v w)  ^ = - (w  r + 2wV  + "u^'r)  + — C F 


(12) 


and,  of  course, 


r = (u  v w)T  (13) 

o 

as  the  equations  which  govern  the  motion  of  m (assuming  that  co  and  u)  are 
known) . 

The  usefulness  of  vector-matrix  conversion  is  now  evident.  General 
digital  computer  subroutines  can  be  written  to  form  tilde  matrices  from 
vectors  and  multiply  matrices  and  hence  evaluate  the  right-hand  side  of 
Equation  (12),  which  will  be  a column  matrix.  Thus,  the  equations  of 
motion  of  a complicated  dynamical  system  do  not  have  to  be  completely 
reduced  to  scalar  form  by  hand  to  program  them. 
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The  advantages  of  using  a digital  computer  to  perform  a great  deal 
of  the  algebra  are  many,  but  two  of  the  foremost  ones  are  that  mistakes 
in  algebra  (especially  minus  signs)  may  be  avoided  and  that  the  programm- 
ing statements  may  be  more  concise,  hence  avoiding  programming  errors. 

In  what  follows,  we  will  use  the  tilde  notation  just  explained  as 
well  as  (1)  an  underline  for  either  a vector  or  a column  matrix,  (2)  a 
double  underline  for  a second-order  tensor  or  a square  matrix,  (3)  a 
superscript  T to  denote  the  transpose  of  a matrix,  and  (4)  a superscript 
-1  to  denote  the  inverse  of  a square  matrix.  Three  other  convenient 
definitions  which  will  be  utilized  are 


and 


10  0 


E 


0 10 


0 0 1 


10  0 


Si 


0 0 0 
0 0 0 


E23  - 


0 0 0 
0 10 
0 0 1 


(14) 


(15) 


(16) 


c . Dynamics 

The  approach  used  in  obtaining  the  equations  of  motion  of  the 
launcher/rocket  system  is  as  follows.  First,  the  vector  equation 
governing  the  rotational  motion  of  the  system  about  0 is  obtained.  When 
the  rocket  separates  from  the  launcher,  certain  terms  may  be  dropped 
from  this  equation  to  obtain  the  governing  equation  for  the  launcher 
alone.  Second,  the  equations  for  rotation  and  translation  of  the  rocket 
are  developed.  Two  forms  of  the  equation  of  rotational  motion  of  the 
rocket  are  developed,  since  during  detent  and  tipoff  it  is  convenient  to 
sum  moments  about  the  point  a on  the  rocket  centerline  (Figure  2) . 

(1)  System  Rotational  Motion.  Consider  Figure  5 in  which  the 

system  of  launcher  plus  rocket  is  depicted  as  two  rigid  bodies.  FigureS 

shows  the  vectors  R,  S , and  u,  which  serve  to  locate  the  centers  of  mass 

-a  - 
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Figure  5.  Geometric  description  of  position  vectors. 


of  the  two  bodies,  and  the  generic  position  vectors  jj  and  to  ele- 
ments of  mass  in  the  launcher  and  rocket,  respectively,  from  their  mass 
centers . 

By  definition,  the  angular  momentum  of  the  system  about  0 is 


, . L 

o J 

M 


eL)  X (R  + £l)  dM  + J (R- + -s* + u- + 

m 


/■ 


X (R  + S 4-  u + n ) dm 
— —a  — — m 


(17) 


The  expression  for  Hq  may  be  simplified  by  using  the  definition  of  the 
center  of  mass  of  a rigid  body.  We  have 

M M 
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1 


and 


/ £ m dm  = /i 


dm  = 0 


(19) 


m 


m 


and,  since  R,  R,  S_  , S_  , u,  and  u do  not  vary  during  the  indicated  inte- 

a a 

gration  over  the  masses,  we  may  simplify  Equation  (17)  to  obtain 
Equation  (20) : 


u 


H = M R X R + In  X p dM  + m(R  + S + u) 
o J — L — L a 

M 


X (R  + S + u)  + ip  xp 
— —a  — / —m  — m 


/, 


X o dm 


(20) 


Let  ju)  and  0 denote  the  angular  velocities  of  the  launcher  and  rocket, 
respectively,  so  that 


Pl  ' - x Pl 


(21) 


and 


p = n x p 

_m  — _m 


(22) 


By  substituting  Equations  (21)  and  (22)  into  Equation  (20),  it  may  be  shown 
that 

H = I_  / • u + m(R  + S -hu)X(R-hS  + u)  + I • fl  , (23) 

where  anc^  are  the  inertia  dyadics  of  the  launcher  with  respect 

to  point  0 and  the  rocket  with  respect  to  its  center  of  mass  C . 

m 

The  time  rate  of  change  of  H , H , is  equal  to  the  external  torque 
about  0,  T . Explicitly, 


H — IT  / • (*)  + (*)  X IT  # • (*)  + m(R  + S + u) 

- =L/o  — =L/o  — ~ —a  — 

L 


x (r  + s + u)  + i *n  + nxi  • n 

— —a  — -m  — — =m  — 


(24) 


m 


15 


where  we  have  assumed  that  and  u are  expressed  in  the  launcher  basis 

and  and  Q in  the  rocket  basis.  Equation  (24)  is  the  basic  equation 

governing  the  rotational  motion  of  the  system  about  0.  We  must,  however, 
consider  the  rotation  and  translation  of  the  rocket  before  we  can  solve 


Equation  (24)  for  w (in  matrix  form),  since  w may  be  present  implicitly 

in  the  third  term^during  the  detent  phase.  Hje  also  point  out  that  in 
Equation  (24)  and  in  certain  equations  which  follow  a small  circle  above 
a vector  and  an  L or  m underneath  it  indicates  the  derivative  of  the 
components  only  of  such  a vector  expressed  in  the  launcher  (L)  or  rocket 
(m)  basis. 


(2)  Rocket  Rotational  and  Translational  Motion.  Considering 
the  rocket  only,  its  equation  of  rotational  motion  may  be  expressed  as 


i •n  + nxi  • a = t , 

-m  “ — -m  — —m/cm 

m 


or  as 


(25) 


I (1  + fi  X I • fl  = -u  X m(R  + S + u)  + T / 

— — —a  — —m/a 


(26) 


=m 


m 


Where  — m/C  and  are  the  torques  on  the  rocket  about  its  center  of 


m 

mass  and  about  a,  respectively. 

The  translation  of  the  rocket* s center  of  mass  is  governed  by  the 
equation, 


m(R  + £>  + u)  = F 

— —a  — — i 


(27) 


where  F is  the  external  force  on  the  rocket.  More  explicitly,  we  have 
— m 


/A  /A  /N  /A 

'‘That  is,  I , = I i i - I i j - I i 

=L/o  x x L L x y L L x z 

Li  Li  L Li  Li  Li 


A A A A /N  A A 

+ I jTjT  - I jTkT  - I kT  i - I k j 
yLyL  L L yLZL  L L XLZL  L L yLZL  L L 


AAA 

U)  = + + simi-lar  definitions 


£ - i 3 £ 
Vl  L L 

H-  I k_  k_  and 

z z L L 
L L 

for  I and  Q. 

=m  — 
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m 


u>XR+u>X(u>XR)+S  + 2 u)  X S + u)  X S 


f L 


+ uix(u)xs)  + nxu  + nx  (axu) 

m 


= F 


"m 


(28) 


Equations  (24)  through  (28)  may  be  used  in  a variety  of  ways.  What 
we  really  want  to  achieve,  of  course,  is  a consistent  set  of  equations 

o o oo 

which  may  be  solved  for  w,  ft,  and  S . 

L m —r 


Using  Equations  (26)  and  (27),  Equation  (24)  can  be  rewritten  in 
the  form, 


IT  / • w + w X 

=L/o  — — 

L 


IT  / * w + 

=L/o  — 


(R  + S 
— —a 


u)  X F + T , = T 

— — m —m/cm  — o 


(29) 


which  is  a suitable  equation  for  U),  if  F and  T y_, 

mm/C 


do  not  include  con- 


m 

straint  forces  and  torques;  whereas,  if  they  do,  £>  will  also  appear  in 
those  forces  and  torques.  We  shall  consider  constraints  on  the  rocket 
and  their  effects  on  the  forms  of  the  equations  of  motion  in 
Subsection  Ill.f. 


Before  considering  the  kinematical  aspects  of  our  problem,  we  note 
that  | w|  should  be  very  small  for  most  launchers  so  that  we  may  drop  the 
term  jx)  X i^/0  • from  Equation  (29)  and  similar  terns  which  appear  in 

Equation  (28) . * 

d.  Kinematics 


Kinematics  is  a subject  area  which  deals  with  the  problem  of 
describing  position,  velocity  and  acceleration  in  terms  of  geometric 
variables.  In  our  case,  we  will  use  kinematical  principles  to  (1)  express 
and  ft  in  terms  of  Eulerian  angles  and  their  derivatives,  (2)  obtain  an 
alternate  equation  of  motion  for  the  rocket1 s center  of  mass,  (3)  derive 
equations  for  the  time  rates  of  change  of  X,  Y,  and  Z,  the  inertial 
coordinates  of  C , (4)  obtain  equations  defining  the  orientation  of  the 

rocket  with  respect  to  the  launcher,  and  (5)  develop  the  transformation 

matrices  needed  to  convert  the  dynamical  (and  some  kinematical)  equations 

to  the  (t,3,£)  (inertial),  (i  t ,1),  or  (l  ,j  ,£  ) basis. 

L L L m m m 


*These  terms  have  been  retained  in  the  computer  code. 
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(1)  Launcher  Orientation  and  Angular  Velocity.  The  orienta- 
tion of  the  launcher  may  be  defined  by  the  set  of  Eulerian  angles  (Q^, 

OL^yOL^)  shown  in  Figure  6.  The  rotation  sequence  of  2-3-1  was  chosen  so 


Y 


L 


Figure  6.  Launcher  orientation. 


that  large  values  of  the  elevation  angle  maY  be  easily  incorporated. 

It  is  assumed,  however,  that  ^ anc*  OL ^ are  small  angles  so  that 

With  this  assumption. 


cos  *2  ~ cos 


1,  sin  ^ ^ a^,  and  sin 


the  direction  cosine  matrix  which  relates  the  (I,J,K)  and  (i  , j ,£^) 


AAA 


bases;  i.e.,  (i  j k )T  = B(I  J K)T  is 

L L L 
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B ~ 


ca„  a„ 


OjScr2  - a3ca2 


-sa„ 


1 a3Sa2  + a1QC0:2 


sol  -a. 


ca„ 


(30) 


where  the  shorthand  notation  S = sin  and  C = cos  has  been  used. 


AAA 

The  components  of  w = w i H~  (aL  j + U)  k , may  be  expressed  in 
. r 1 L Z L j L 

the  forms, 


= a - a2  s <3 


(31a) 


w2  = a2  ca  co^  + a3s  l 


(3lb) 


and 


W3  = "a2CO,3SO!l  + a3Cai 


(31c) 


Since  the  CC.  are  required  to  be  small  and  a and  Q t are  small  also, 


tx)  ~ a 

1 1 


(32a) 


and 


U)A  ^ 


ql 


y 


or 


(32b) 


(32c) 


(33) 


(2)  Rocket  Orientation  and  Angular  Velocity.  Because  the 
rocket  is  assumed  to  be  connected  to  the  launcher  by  springs,  the 
orientation  of  the  rocket  relative  to  the  launcher  as  well  as  its 
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inertial  orientation  should  be  known.  We  also  require  knowledge  of  the 
rocket's  angular  velocity  relative  to  the  launcher,  since  the  springs 
possess  structural  damping  characteristics. 

The  inertial  attitude  of  the  rocket  may  be  defined  by  a 3-2-1 
sequence  of  Eulerian  angles  denoted  by  \|/,  0,  and  <t>  , respectively,  and 
displayed  in  Figure  7.  The  matrix  Q which  defines  the  transformation 


Figure  7. 


Inertial  orientation  of  the  rocket. 


AAA  AAA 

from  the  (I  J K)  basis  to  the  (i  j k)  basis  is,  in  rotation  matrix 


‘i 

0. 

0 ‘ 

’ce 

6 

-se" 

C\|r 

sy 

o' 

Q = 

0 

c<t> 

S4> 

0 

1 

0 

-S\|/ 

0 

0 

-S<t> 

c<t> 

_se 

0 

ce  _ 

_0 

0 

1 

AAA 


and  the  components  of  Q = 
forms , 


H,  i + (]„  j + k may  be  expressed  in 
lm  2 Jm  3 m r 

^ M - t S G 

d2  = 0 c$  + f ce  s<t> 


form. 


(34) 


the 


(35a) 

(35b) 


and 

n3  = -e  s<t  + ijr  ce  c<i>  . (35c) 

Equation  (35)  may  be  used  to  obtain 


0 = Q.l  + (ft2S<t>  + fl  C«)  tan  6 


(36a) 


e = a2  c<t>  - n3s  <t> 


(36b) 


t = 


n2s<t>  + ft3  c<t> 
ce 


(36c) 


Since  only  the  angle  may  be  assumed  to  be  small.  Equation  (36)  may  not 
be  simplified. 

The  rocket's  orientation  relative  to  the  launcher- fixed  C x y z 

L L L L 

system  may  be  defined  by  the  Eulerian  angle  set  (0^,0O,0^).  These  angles 

are  depicted  in  Figure  8.  The  rotation  sequence  is  the  same  as  that  for 
the  set  (\|/,0,<t>)  but  both  0^  and  0 ^ may  be  assumed  to  be  small.  By  using 

the  angles  0 we  may  form  a matrix  A which  specifies  the  transformation 

from  the  (iT,jT,tt_)  basis  to  the  (i  j k ) basis.  Explicitly, 

L L L m m m 
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Figure  8.  Orientation  of  the  rocket  relative  to  the  launcher. 


-6, 


A ^ 


+ ^S^i  C 6^ 


63SO1  + Ce 


(37) 
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The  angular  rates  6 are  related  to  Q and  u)  by  the  equation, 


- e3  se2 


Q - A u>  = 


e2cel  + e3ce2se1 


-e2 sel  + e3  ce2  cel 


Since  Q 2 and  6 3 are  small. 


1 

O 

o 

1 

“1 

A 'a)  - 

o ce,  se. 

U)_ 

l i 

2 

o -se,  ce. 

UK 

l l 

L 3 

so  that 


^1  ~ ^1  + 


o)3)tan  02  - w 


and 


e2  « a2  ce  - a3  se  - w2  , 


S3  ~ 


q2  S01  + n3  ce  - oo3 
cel 


(38) 


(39) 


(40a) 

(40b) 

(40c) 


Assuming  that  0 „ is  appropriately  small, 


9, 

fi,  - U), 

1 

1 1 

0O 

CZ 

n„ce,  - o„se,  - 

2 

2 1 3 1 2 

fi2se.  + n ce^^  - w3 

L 3^ 

(41) 


Equations  (30),  (33),  (34),  (36),  (37),  and  (41)  are  sufficient  to 
connect  the  rotational  motions  of  the  system  with  the  physical  world 
through  geometry.  We  still  need  analogous  equations  for  the  translational 
motion  of  the  rocket. 
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(42) 


Let  the  velocity  V of  be  expressed  as 


V = U i + V i + W k 
— m m m 


We  then  may  write 


Since 


V = V + H X V 


V = 


0 — m 

v = a x v + — 

— — “ m 

A matrix  form  of  Equation  (45)  is 

U 
V 

. AAA 

if  we  assume  that  F is  expressed  in  the  (iT,j_,kT)  basis. 

— m L L L 

The  velocity  of  C may  also  be  expressed  in  the  form, 
m 

• A • A • A 

V=XI  + YJ  + ZK 


= -n  v + a — 

m 


and  it  then  follows  from  the  definitions  of  A,  B,  and  C that 


X 

U 

u 

T T 

T 

Y 

PQII 

II 

V 

= c 

V 

_Z_ 

_w. 

_w_ 

Because  we  also  want  to  know  the  vector  S , we  let 

—a 

' XL  \ + Y1  5L  + ZL  \ • 

a a a 


(A3) 


(AA) 

(A5) 


(A6) 


(A7) 


(A8) 


(A9) 
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5 


(50a) 


XL  =XL  +&1 
a a 

o 


YL  =YL  +&2  ’ 

a a 

o 


(50b) 


and 


ZL  = ZL  + &3  ' 

a a 

o 


(50c) 


where  X , Y , Z are  initial  values  of  the  indicated  variables,  and 

Li  Li  Li 

a a a 
o o o 

consider  Equation  (28).  We  thereby  obtain  the  matrix  equation, 


_ 3_ 


P <rv< 

= — - 5)  (R  + s ) 
m L a 


rsj  / • « • \ 

- V l°2t'3>) 


rp  rp  q rp  iv 

- A ft  u - A ft  ft  u 

m 


(51) 


We  now  have  dynamical  and  kinematical  equations  for  the  system  when 
there  are  no  constraints,  since  Equations  (25)  and  (29)  are  then  uncoupled 
as  far  as  amd  $ are  concerned . Thus,  Equations  (25)  and  (29)  can  be 
solved  for  the  matrix  equivalents  of  these  vectors  and  the  results  sub- 
stituted into  Equation  (51)  when  there  is  no  detent  constraint.  However, 
when  the  rocket  is  constrained  by  the  detent  force,  which  we  assume  acts 
through  the  point  a,  Equations  (25)  and  (29)  are,  in  general,  coupled, 
and  equations  which  account  for  this  fact  must  be  obtained  for  that 
phase  of  the  launch.  We  also  need  to  convert  the  remaining  equations 
for  rotational  motion  into  appropriate  matrix  equations  and  to  obtain 
expressions  for  the  forces  and  moments.  These  things  are  done  in  the 
next  two  subsections. 

e . Constraints  and  Matrix  Forms  of  the  Equations 

The  detent  force  is  assumed  to  be  the  only  nhardM  constraint 
on  the  rocket's  motion;  i.e.,  it  is  the  only  constraint  that  completely 
prohibits  motion.  This  force,  which,  of  course, varies  in  magnitude  and 
direction  during  spinup  and  thrust  buildup,  is,  as  stated,  assumed  to 
act  through  the  point  a,  and  to  prohibit  motion  of  the  rocket  in  the 
x^  direction.  This  constraint  may  be  defined  mathematically  as 

xL  = o (52) 

a 
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and 


xL  = 0 (53) 

a 


during  detent.  If  we  let 


F = F + F + F 
— m — m 

ex 


_mD  ““F 


where  F is  the  resultant  of  the  external  forces  on  the  rocket  (mclud- 

— m 

ex 

ing  those  due  to  gravity,  thrust  and  the  rocket  springs),  F^  is  the 
detent  force  and  F^  is  the  frictional  force,  then  Equations  (52)  and 
(53)  are  satisfied  if 


i 


F -2m 
— m 

ex 


— X (^2  + ^3  X 0}  x 3J) 


+ mu  X H + m(R  + S ) X uj  ~ F +F 

_ _ "a  l _mD 


(54) 


The  detent  force  has  only  an  x -component  in  our  model,  thus 

Li 


_mD 


2 m u X ^ + m £,  X (£  X u) 


- m(R  + S)Xu)-muXfi-F  -F 

1 — o — — — * — m — 


m 

m ex 


•J  ■ '•■! 5 


L * (55) 


The  rocket  will  start  to  move  in  the  x-  -direction  when 

Li 


F 

— m 


ex 


2 m w X (&2  + 63  k)  - m ft  X (ft  X u) 


+ m(R  + S ) X w+muXft 

n_  —a  — — — 

L m 


> 

F 

+ 

F 

_mF 

max 

(56) 


Inequality  (56)  may  be  used  as  a check  to  determine  when  detent  release 
occurs . 
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Consider  now  Equations  (26),  (27),  (54) 
equations,  we  obtain 


and  (55) . Combining  these 


° r ° - 1 

' !l  " u X ™(R  + j3  ) X w ' i 

m L a L LJ 

- u X m u X ft  • i ll, 

m -* 


T + ft  x i - a 

L — =m  — 


- X [(\X  • \)\  + (\x  • \)\] 

!m[-J  X (l2  \ + g3  El) 


u X ( 2 


+ m[n  X (fl  X u)]  • j_  \ i_  + T , 

— x-  — L ( L —m/a  , 


(57) 


as  the  equation  of  rotational  motion  of  the  rocket  before  detent  release. 
In  a similar  manner,  we  obtain  for  the  launcher,  the  equation, 


V>  — ® + ^a)  x ”L<s  + V 

L 


j«,[(R  + Sa)  xi]-  iLji 


L - ® + V 


/ 0 \ ^ 

l„  1 

1 1 

X 

ml  u X n 1 • iLJ 
- N m ' J 

1L  - -®  + sa)  x| 

iH 

[a  x (62  jL  + 63  V 

• 1L  + m[n  X (B  X U)]  • \J  tL  - (2  + sa)  x F(f 

+ /f  • U k,  + T - T , 

I -m  Li  1 -o  —m/a 


ex 


A \ A 
)j> 


JL  Pi 


(58) 


Equations  (57)  and  (58)  may  be  combined  to  form  the  matrix 
equation, 


hi 


hi 


hi 


=22 


-u  A EOQ  F + T / - ft  I 

— = =23  ~m  —m/a  — =m  — 
ex 


(R  + S ) Eoa  F - A1  T , + T 

“a  =23  — m = — m/a  ~ o 


ex 


■2  m u A E go  E „ S_  - m u A E AT  Q S u 
- - -l-  -23  - = =1  - 


-2  m (R  + Sa)  Ex  w g23  Sa 

L 
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m (R  + S ) E.  A ft  ft  u 
a i 


* (59) 


AAA 


where  we  have  assumed  that  T , is  expressed  in  the  (i  ,j  ,k  ) basis  and 

“m/a  m m m 

AAA 

T and  F in  the  (i_,jT,k  ) basis  and  where 
“O  m L L L 


J. . = I - m u A E,  A u 
=11  =m  — = =1  = — 


=12  = "m  5 A it  (5  + , 


J01  s -m  (£~+~S  ) E.  AT  u = 

=21  —a  =1  = “ =12 


and 


=22  s h/o  - “ ® + Sa)  h <S  + Sa> 

Hence,  by  inverting  the  matrix, 

J = 


ill 

1 

1 

=12 

^21 

1 

1 

i22 

and  premultiplying  Equation  (59)  by  j's  inverse’,  L,  we  get 


O 

L, , T + L,  „ T, 

=11  -m  =12  -L 

o 

UJ 

L„,  T + Tt 

=21  -m  =22  -L 

where  T and  T are  defined  such  that 
— m — L 


O 

fi 

T 

m 



“m 

O 

UJ 

L 

(60a) 

(60b) 

(60c) 

(60d) 


(61) 


(62) 


(63) 
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The  inverse  of  J may  be  determined  by  inverting  only  3X3  matrices, 

which  can  be  easily  done  explicitly.  The  formulas  for  computing  the 

L . . are  [6] 

=ij 


=11  (=11  " =12  -22  =21^ 
=22  = (=22  ' =21  =11  ^12) 


(64a) 


(64b) 


=12  "=11  =12  =22 


(64c) 


and 


=21  "-22  =21  =11 


(64d) 


The  matrix  equations  of  motion  for  post-detent  phases  are 


and 


I 

=m 


(]--(]  I 1]  + T , 

— — =m  — m/Cm 

m 


T 

AT, 

= — m/Cm 


where 


(65a) 


(65b) 


F = F + F 

-m  -mex  -nip 

Equations  (30),  (33),  (36),  (37),  (41),  (46),  (48),  (51),  and  (62) 

or  (65)  are  those  required  to  simulate  the  motion  of  the  system.  However, 

F , T / .and  T must  be  defined  before  these  equations  may  be  solved. 

— m “m/a  ~ o 

Such  definitions  are  given  in  the  following  subsection. 

f . Forces,  Torques,  and  Equilibrium 

(1)  Forces  and  Torques  on  the  Rocket.  The  forces  which  act  on 

the  rocket  were  assumed,  for  the  purposes  of  this  study,  to  be  (1)  the 

force  due  to  gravity,  F (a  flat  earth  is  assumed),  (2)  the  thrust  force, 

m 

g 

F^  , (3)  the  detent  force,  F^,  which  has  previously  been  discussed, 

(4)  a frictional  force  which  opposes  the  rocket1 s motion  parallel  to  the 
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x -axis,  and  (5)  the  forces  F and  F in  the  fore  and  aft  springs  which 

L IS  HS 

couple  the  rocket  to  the  launcher.  Explicit  expressions  for  these  forces 
are  derived  in  what  follows. 


(a)  Gravity  Force. 
mgK  so  that  in  matrix  form  (in  the 


The  gravity  force  on  the  rocket  is 

iL.3L»^L  basis) 


0 


F = B 0 
— m = 

g 

mg 


(66) 


where  g is  of  course  the  local  acceleration  of  gravity. 


(b)  Thrust  Force. 


The  thrust  force  F is 

"“t 


obtained  by 


assuming  that  the  thrust  is  rocket- fixed  and  is  misaligned  with  the 

x -axis  only  a small  amount.  As  shown  in  Figure  9,  the  angular  thrust 
m 

malalignment  angles  are  denoted  by  OL  (a  rotation  about  a line  parallel 

to  the  y -axis)  and  a (a  rotation  "essentially"  about  a line  parallel 
m z 

to  the  z -axis, 
m 


Figure  9.  Thrust  vector  definition. 


The  magnitude  of  _F  is  denoted  F,  so  that, 


F_.  ^ F i + Ct  Fj  - Of  F k 
— T m z m y m 


(67) 


From  Equation  (67),  it  follows  that,  in  matrix  form, 


f = A .Fryi 
■ T 


(68) 
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(c)  Spring  Forces.  The  forces  in  the  springs  which  couple 

the  rocket  to  the  launcher  are  assumed  to  be  linear  functions  of  the 

relative  displacements  and  velocities  in  planes  parallel  to  the  y z -plane 

L L 

of  the  points  f and  a,  measured,  of  course,  from  their  undeformed  posi- 
tions in  these  planes.  Thus,  the  force  in  the  aft  springs  is  (matrix 
form) 


"0 

“0  " 

F = -C 

-K 

—as  =a 

2 

a 

2 

-&3_ 

(69) 


where  C and  K are  damping  and  stiffness  matrices,  respectively, 
a a 

To  determine  the  force  in  the  forward  springs,  we  consider  Figure  10. 
From  this  figure,  we  conclude  that 


"0  “ 

‘0  “ 

*2 

+ A 

>3. 

A 

"0  " 

"0  ' 

&2 

+ ifa 

CM 

CD 

L 3 J 

L 3 J 

(70) 


Figure  10.  Displacements  of  points  f and  a. 
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where , 


0 0 0 

0 0 -l 

0 i,  0 

fa 


(71) 


and  £ is  a negative  number,  since  it  is  the  x -component  of  £c  which 
ta  m — la 

lies  on  the  rocket  centerline  and  is  directed  from  f to  a. 

(d)  Frictional  Force,  Guidance  mechanisms,  such  as  rails 

and  shoes,  inherently  are  friction  producers.  The  frictional  force  on 

the  rocket  opposes  motion  and  in  our  model  is  assumed  to  act  only  parallel 

to  the  x-axis.  Since  the  spring  forces,  F-  and  F , act  perpendicular 
Li  t S 3.S 

to  the  x^-axis,  it  is  assumed  that  the  frictional  force* s magnitude  is 

p(  | F | + |Ff  |),  where  p is  a friction  coefficient.  The  direction  of 
a s t s 

the  frictional  force  is  assumed  to  be  in  the  -x^-component  of  the 
rocket's  weight.  Therefore, 

+ llj)lL  . (72) 


(e)  Torque  About  C^.  The  torque  about  is  due  to  the 

thrust,  the  springs,  and  friction.  Explicitly,  in  vector  form,  we  have 
for  center-shoe- type  rocket  guidance, 


= - (u  - Av)  X IT  - (U  + i,  ) x F. 


- u xIas  - (-“  - Afa±  Pf  jL)  X ^ Ilfs  IS, 


(-U  ± Pa  JL)  x p|l,J  iT  + (!  - p0)l 


as  L 


s — s 


(73) 


where  and  are  the  distances  to  fore  and  aft  rails,  respectively, 

from  the  rocket  centerline  and  the  signs  on  p and  p are  determined  by 

t a 

the  signs  of  the  y -components  of  Ff  and  F so  that  the  frictional 
Li  r s as 

forces  act  on  either  the  right  (+)  or  left  (-)  side  of  the  rocket  at  the 

fore  and  aft  shoe  positions.  Also,  in  Equation  (73),  £ , a vector  from 

Jr 

a to  a point  on  the  line  of  action  F may  be  expressed  as, 
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(74) 


= X i -fa  j + a k 

— F F m y m z 


where  a and  a are  normally  much  less  than  the  rocket !s  diameter  and 

y z 

represent  the  amount  of  linear  thrust  malalignment.  Finally,  \i  is  the 
spin-torque  friction  coefficient  and  is  the  spin  torque. 

A suitable  matrix  form  for  T (rocket  basis)  is 

“m/C 

m 


T , = -(u  - i)  F ■ (U  + I ) A F.  - 

-m/C  — F — T -fa  = -fs 

m 

+ (U  + jLfa  " Pf)  P llfg  l d00)T 

+ Cm  ■ Pa)  p|Fasl(i°o)T  + (l  - ns) 


u A F 
— = —as 


T 

— s 


> 


(75) 


where 


pf  = 1 pf 


0 0 1 

0 0 0 

-1  0 0 

0 0 1 

0 0 0 

-1  0 0 


(76) 


(77) 


Tn/a’ 


(f)  Torque  About  a.  The  torque  on  the  rocket  about  point  a, 
may  be  simply  expressed  in  the  matrix  form, 


T / 
—m/a 


Tl/C 


m 


u A F 
— - — m 


(78) 


(2)  Torque  on  the  System.  The  torque  is  that  about  0 due 

to  external  forces  on  the  system  of  launcher  plus  rocket.  We  assume 
that  the  only  external  forces  acting  on  the  system  are  those  due  to 
gravity,  the  rocket's  propulsion  system  and  the  launcher  supporting 
structure  which  is  modeled  as  pitch,  yaw,  and  roll  torsional  springs. 
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It  is  assumed  that  the  natural  frequencies  of  the  launcher's  vibra- 
tion in  pitch,  yaw,  and  roll  are  known.  We  designate  these  natural  fre- 

, respectively,  and  the  associated 


quencies  (undamped)  by  u)  , u) 


n n 

yT  z 


CO 


damping  ratios  by  £ , t , 

yL  ZL 


L “L 
and 


We  also  assume  that  the  launcher 


equations  of  motion  can  be  written  as 


iT  a + cT  a + kt  a = f , 

=L  o — =L  — =L  — — 


(79) 


where  _f  may  be  a nonlinear  function  of  Ct  = (CH^  CX ^)  and  its  deriva- 

tives and  any  other  of  the  variables  which  define  the  system  motion,  and 


also  as 


so  that, 


a = 2 


L 0 0 


o r o 

yL 


0 0 


U) 


XT 


U) 


U) 


a 


+ 


U) 


U) 


U) 


a = i_  / f 

— =L/o  — 


(80) 


CT  = IT  , 

=L  =L/o 


2t  u) 

"xT  n 
L x. 


it  u)  o 
yT  n 

L yL 


It  U) 

^zT  n 
L z. 


(81) 
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and 


£l 


Il/o 


0) 


XT 


0 0 


U) 


0 0 


0) 


(82) 


are  the  launcher  damping  and  stiffness  matrices. 

The  torque  about  0 may  now  be  expressed  in  the  matrix  form  (launcher 
basis) , 


T = -CT  a - Kt  a + ( R + S AT  u A ) B mg  + R B Mg 
— o *L  — =L  — \—  — a « — =/  = = “ * = 

~ T ~ \ T T 

R + S + A £ A ) A F + A T 
— —a  = — F = / = — T = — s 


(83) 


(3)  Equilibrium  Considerations,  It  is  assumed  that  initially 
the  system  is  in  equilibrium  under  gravity  and  elastic  structural  forces. 
It  is  further  assumed  that  the  small  elastic  displacement  of  the  rocket 
relative  to  the  launcher  may  be  neglected  in  determining  the  equilibrium 
orientation  of  the  launcher,  but  that  the  equilibrium  orientation  of  the 
launcher  should  be  used  in  computing  the  rocket1 s relative  equilibrium 
displacement . 


(a)  Launcher  Equilibrium.  Let  Q denote  the  orientation 

— 0 

vector  of  the  launcher  when  the  launcher  springs  are  not  deformed  and 
let  Of  denote  the  change  in  a which  will  result  in  a balance  in  the 

spring  and  gravitational  forces.  The  equilibrium  condition  then  is 

“(I  + Ia  + ii>!e  & + M 1 ^ ae  , (84) 

where,  without  loss  of  generality,  it  is  assumed  that  A = E initially 
and  where  for  small  \cx 

e i 


B s*  -3  B + B 
=e  — e -o  =o 


(85) 


with  g denoting  the  matrix  g,  in  the  undeformed  state. 
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By  using  Equation  (85)  in  Equation  (84)  > we  obtain 

j[m(R  + la  + u)  + MR]  (jT"g)  - K^j  = Z , (86) 

where 

Z = - (m(R  + 3 + u)  + MR]  B g . (87) 

— ^ —a  — — =o  — 

Therefore, 

-e  = '{&L  ' [m<^  + + ^ ‘ ^]  (D>}_1  % (88) 

and  the  initial  value  of.  Ct  - CL  + OL  . 

— — o “ e 

(b)  Rocket  Equilibrium.  Since  the  rocket  may  translate 
in  the  y - and  z -directions  (x  —motion  constrained)  as  well  as  rotate, 

Li  Li  Li 

two  vector  equilibrium  conditions  must  be  satisfied;  i.e.,  the  sum  of 
the  forces  acting  on  the  rocket  must  be  zero  and  the  sum  of  the  torques 
about  a point  on  the  rocket,  say  point  a,  must  also  vanish. 

Summing  forces,  we  have 

Fr  + F + m B g = 0 (89) 

— fs  —as  =e  — — 

and  summing  torques  about  a,  we  find  that 

-2fa  Ifs  + m u Be  g = 0 . (90) 

It  follows,  from  Equations  (89)  and  (90),  that 

—fa  -as  = ‘m(^  + — fa^=e  £ * (91) 


Let 


mgb  = m(u  + Jfa)|e  g 
Then,  the  equilibrium  values  of  an8  are 


-mgb. 


2 x , k 
e fa  a 


33 


(92) 


(93) 
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and 


mgb 


3 x,  k 
e fa  a 


(94) 


22 


where  xr  i = £-  and  k and  k are  the  elements  of  K indicated  by 
fa  m —fa  a^2  =a 

the  subscripts. 

Once  5^  and  &3  are  known,  6 and  6^  may  be  obtained  from 


and 


-cr 


x£a  kf 


33 


f 


22 


(95) 


(96) 


where  c ^ and  are  the  indicated  elements  of 


£ = ” if,  5e  s - —fa  Sf 


(97) 


g , Rocket  Inertia  Matrix 


The  centroidal  inertia  matrix  of  the  rocket,  L , is  obtained 

=m 

by  assuming  that  the  principal  centroidal  inertia  matrix,  I , is  known. 
Let  ^m 


(98) 
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and  let 


1 


D = 


‘^3 


ll2 


^3 

1 

0 


0 

1 


(99) 


denote  the  matrix  required  to  transform  from  the  C x y z system  to  the 

m m m m 

principal  system  C^x^y^z^  (Figure  11) . Then,  1^  is  given  by 


T 

I = D I D 

=m  = =p  = 

m 


(100) 


Figure  11.  Rocket-fixed  and  principal  coordinate  systems. 


38 


which,  neglecting  squares  and  products  of  ^ an<^  ^3*  explicitly  is 


I - 
=m 


^3^1  ’ 12)  ’ *3^ 


^3a1  - I2) 


- 2 ( I i " I3') 


| Since  we  may  express 

I 

=m 

as 

I 

XX 

-I 

xy 

-I 

xz 

I 

=m 

= 

-I 

xy 

I 

yy 

-I 

yz 

-I 

xz 

-I 

yz 

I 

zz 

it  follows  from  Equations  (101)  and  (102)  that 


1 ~ It  » 

XX  1 

*xy  '-*‘3<I2  ' V 

J 

1 * l0  > 
yy  2 

“ •*2<I3  ' V 

> 

I ^ L and  I « 0 
zz  3 yz 


(101) 


(102) 


39 


Section  IV.  SIMULATION  METHODS 


a.  General  Comments 


The  equations  which  govern  the  motion  of  the  system  have  been 
completely  defined  except  for  the  magnitudes  of  the  spin  torque  and 
thrust.  Since  the  manner  in  which  the  time  variations  of  these  quanti- 
ties is  generated  is  more  of  simulation  topic,  it  is  therefore  treated 
in  this  section.  This  section  also  discusses  matrix  manipulations  used 
in  generating  the  right-hand  sides  of  the  time  derivatives  of  the  state 
variables  and  the  method  of  integration  used  to  solve  the  set  of  27 
differential  equations. 

b . Generation  of  Spin  Torque  and  Thrust 

The  magnitudes  of  the  spin  torque  and  the  thrust  are  generated 
for  computational  purposes  in  the  same  manner;  hence,  we  will  describe 
the  procedure  for  the  thrust  magnitude,  F. 


Values  of  F at  times  when  derivatives  are  needed  are  obtained  by 
assuming  that  the  thrust  profile  is  piece-wise  linear.  In  the  present 
version  of  the  simulation  program,  four  values  of  the  thrust  magnitude, 
Fy  j = 1,2, 3, 4,  at  times  t^,  j = 1,2, 3, 4,  are  used  for  this  purpose 

(Figure  12).  When  T^  - t < t^^,  the  value  of  F is  found  from 


F - F.  + (F  - - F . ) (t  - t.Ht.,.  - t.) 
J j+1  J J j+1  j 


(103) 


For  t < t^  (in  the  computer  code,  T5  = t^,  T6  = t^,  ...  T9  = t,_)  , F = 0 
and  for  t > t^ , F - 0. 

c . Matrix  Algebra 

The  matrix  operations  which  are  needed  to  generate  the  deriva- 
tives used  in  the  integration  process  are  (1)  multiplication  of  a 3 X 3 
matrix  times  a 3 X 1 matrix,  (2)  multiplication  of  two  3x3  matrices, 

(3)  generation  of  a "tilde"  matrix  from  a 3 X 1 matrix,  (4)  inversion 
of  a 3 X 3 matrix,  and  (5)  matrix  addition,  subtraction  and  transposition 
Matrix  subtraction,  addition,  and  transposition  are  carried  out  directly 
when  necessary  during  the  computations  by  using  DO  loops.  The  other 
operations,  or  manipulations,  are  carried  out  in  separate  subroutines, 
all  of  which  are  very  straightforward  because  all  elements  of  the  result 
are  calculated  explicitly  even  in  the  matrix  inversion.  For  completeness 
the  algorithms  of  the  subroutines  MATXV  (3x3  times  3x1),  MATMPY 
(3x3  times  3x3),  TILDE  (form  3x3  matrix  from  3x1),  and  MATINV 
(inverse  of  3 X 3)  are  given  as  follows : 1 


tThese  operations  are,  of  course,  well  known. 
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THRUST  MAGNITUDE,  F 


TIME,  t 

Figure  12.  Thrust  magnitude  versus  time. 


MATXV 


Input: 


11 

a12 

a13 

21 

a22 

a23 

• 

(104) 

31 

a32 

a33 

Operations: 


W . = L a k vk,  k = 1,2,3,  j = 1,2,3 
J k J 


(105) 


Output: 


w 


1 


w 


w 


3 


(106) 
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MATMPY 


\ 


Input : 


Operations : 


Output : 


TILDE 


Input : 


Output: 


all 

a12 

a13 

4 = 

a21 

a22 

a23 

> 

a32 

a33 

b12 

b13 

B = 

b21 

b22 

b23 

b31 

b32 

b33 

(107) 


(108) 


c 


ij 


I 


k=l 


ik 


(109) 


11 

C12 

C13 

21 

C22 

C23 

• 

(HO) 

31 

C32 

C33 

X = 


X1  X2 


(111) 


X = 


(112) 
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MATINV 


Input: 


Operations : 


all 

a12 

a13 

IO 

II 

a21 

a22 

a23 

a31 

a32 

a33 

all 

(a22 

a33  “ a32  a: 

a12 

(a31 

a23  " a2i  a: 

a13 

(a21 

a32  " a31  a 

11 


<a22  a33  ’ a32  a23^° 


12 


^32  a13  ’ a12  a33^/D 


13 


(a12  a23  ’ a22  a13)/D 


21 


^a31  a23  ’ a21  a33^D 


22 


(all  a33  “ a31  a13)/D 


23 


(a21  a13  ‘ all  a23)/D 


31 


^a21  a32  " a22  a31')  J 


32 


^a31  a12  " a32  all)/D 


33 


(all  a22  ' a12  a21)  ° 


Output : 


B , DET  = D 


Inversion  of  the  6x6  matrix  J has  already  been  discussed  in 
Subsection  Ill.f. 


(113) 


(114) 


(115) 
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d. 


Integration  Method 


Various  methods  may  be  used  to  solve  sets  of  first-order 
ordinary  differential  equations.  The  most  popular  type  of  self-starting 
algorithm  is,  however,  the  Runge-Kutta,  and  a fourth-order  algorithm  of 
this  type  was  chosen  for  use  in  this  study. 

The  fourth-order  algorithm  is  basically  as  follows  [7]: 


Given  the  state  vector  x (a  27  X 1 matrix  in  this  case)  at  time  t 
and  the  time  derivative  functions  f ^ (x,t);  i.e., 


x = f (x,t)  , 

1)  Compute: 

wx  = f(x(t),t)At 

where  At  is  the  step  size. 

2)  Compute: 

u^  = x(t)  + w^/ 2 

3)  Compute: 

w^  = f(u^,t  + At/2)At 

4)  Compute: 

u 2 = x(fc)  + 2J2/2 

5)  Compute: 

= J(u2»t  + At/2)At 

6)  Compute: 

u^  = x ( t ) 4-  w^  . (122) 

7)  Compute: 

w4  = f(u3,t  + At)  At  . (123) 

8)  Compute  the  new  state  vector,  x(t  + At) , from 

x(t  + At)  = x(t)  + + 2— 3 + • (124) 

By  redefining  t to  be  t + At  and  x(t)  to  be  x(t  + At)  , one  may  repeat  steps  1) 
through  8)  until  t is  as  large  as  desired , or  until  one  of  the  elements  of  x has 
attained  a certain  value. 


(116) 

(117) 

(118) 

(119) 

(120) 
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Section  V.  SIMULATION  RESULTS 


a . General  Comments 


As  a check  on  the  theory  and  the  computer,  several  test  runs 
were  made  using  data  which,  it  is  emphasized,  do  not  represent  any  actual 
launcher/rocket  system,  but  some  of  which  does  correspond  to  physical 
data  for  the  AERO  Wind  Sensitivity  Rocket  currently  awaiting  tests  by 
the  US  Army  Missile  Command  [8].  The  data  used  to  obtain  the  results 
presented  here  are  given  in  the  appendix  as  an  example  of  the  input 
needed  for  the  computer  program. 

Only  a few  of  the  more  interesting  plots  which  were  obtained  during 
this  study  are  presented  in  the  next  subsection. 

b . Discussion  of  Some  Typical  Results 

Plots  of  various  variables  are  shown  in  Figures  13  through  28. 
The  plots  were  made  and  are  exhibited  to  illustrate  the  capabilities  of 
the  program  and  also  some  of  the  effects  of  the  mallaunch  factors: 

(1)  detent  force,  (2)  mass  unbalance,  and  (3)  thrust  malalignment. 

Five  simulation  runs  were  made.  These  were  obtained  in  the  follow- 
ing manner: 

1)  Run  1 — No  detent  force,  mass  unbalance,  or  thrust  malalignment 
were  included. 

2)  Run  2 — A detent  force  of  1800  pounds  (maximum  magnitude)  was 
applied . 

3)  Run  3 — The  detent  force  of  Run  2 was  retained  and  a mass 
unbalance  caused  by  rotating  the  principal  axis  of  the  rocket  through 
anangle  ^ = 0.00001*  radian  (about  the  y^-axis)  was  included. 

4)  Run  4 — The  detent  force  was  retained  and  a thrust  malalign- 
ment due  to  a = 0.0001x  radian  was  included. 

y 

5)  Run  5 —The  detent  force,  mass  unbalance  of  Run  3,  and  thrust 
malalignment  of  Run  4 were  all  included. 

Figures  13,  14,  and  15  show  the  time  histories  of  the  launcher  roll, 
pitch,  and  yaw  angles,  respectively.  Since  the  launcher  is  relatively 
stiff  and  massive  (compared  to  other  possible  launchers),  the  disturbing 
influences  of  mass  unbalance  and  thrust  malalignment  on  the  launcher 
motion  are  negligible  compared  to  the  effects  of  the  detent  force. 


^This  is  relatively  small  compared  to  values  which  may  be  expected. 
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Note  that  holding  the  rocket  during  spin-up  caused  a noticeable  varia- 
tion in  launcher  roll  and  yaw.  However,  noticeable  is  probably  a 

-3 

misnomer  if  angles  of  less  than  10  radians  are  considered  negligible. 
The  detent  force  has  a noticeable  effect  on  roll,  pitch,  and  yaw,  being 
sufficient  to  cause  a pitch  angle  difference  of  about  0.00005  radian 
(amplitude) . 

Figures  16  and  17  show  the  body- fixed  (missile- fixed)  angular 
velocity  components  and  0^  as  functions  of  time.  The  thrust  comes 

on  when  4>  — 3.66  radians  and  the  missile  leaves  the  launcher  when 
<t>  — 13.25  radians,  and  abrupt  changes  in  the  slopes  of  an<^  ^3  are 

noticeable  at  these  values  of  <t>  . The  choice  of  mass  unbalance  and  thrust 
malalignment  happened  to  be  such  that  they  tend  to  cancel  each  other 
until  the  thrust  cuts  off  when  <t>  - 28  radians.  Even  after  thrust  cutoff, 
the  mass  unbalance  effects  persist,  but  for  Run  5 the  oscillations  in 
^2  and  are  of  smaller  amplitude  than  their  Run  3 (mass  unbalance  only) 

counterparts . 

The  variables  and  are  plotted  in  Figure  18.  This  figure  is 

a little  cluttered,  but  the  almost  identical  traces  for  the  last  four 
runs  until  the  rocket  leaves  the  constraint  of  the  launcher  can  be  seen. 
After  leaving  the  launcher  in  Run  3,  a circle"  centered  at  (-0.00105,  0) 
is  described,  indicating  a bias  in  angular  rate  °f  -0.00105  rad/sec. 

During  Run  4,  smaller  circles  located  approximately  at  (0.0008,  0)  and 
(0,  0)  are  described.  The  combined  effects  of  mass  unbalance  and  thrust 
unbalance  result  in  the  curve  labeled  5 which  ultimately  (after  thrust 
cutoff)  produces  the  small  circle  with  center  at  (-0.00105,  0). 

Figures  19  and  20  illustrate  the  behavior  of  the  missile  pitch  rate 
and  yaw  rate  for  the  five  cases.  The  oscillations  in  6 and  for  Runs  3 
and  4 are  clearly  180°  out  of  phase  (as  they  should  be)  and  for  Run  5, 
the  effects  of  mass  unbalance  and  thrust  malalignment  almost  cancel  out 
until  thrust  cutoff. 

Pitch  rate  is  plotted  versus  yaw  rate  in  Figure  21.  The  curves 

shown  are  similar  to  those  in  Figure  18,  but  the  ellipses  described 

after  launch  by  the  curves  for  Runs  3,  4,  and  5 are  observed  to  wobble; 

i.e.,  the  curves  do  not tdoseas  those  of  Figure  18.  This  is  due  to 

long-period  changes  in  0 and  resulting  from  the  small  difference  in 

the  spin  rate  and  the  rate  of  precession.  This  difference  divided  by 

the  spin  rate  is  approximately  I /I  (assuming  I = I ). 

x x y y v yy  zz 

mm  mm  y mm  mm 


*N o t e the  scale  difference  in  f and  ^ which  result  in  circles  rather 
than  ellipses  on  Figure  18. 
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Curves  showing  rocket  pitch  angle  versus  rocket  yaw  angle  are  pre- 
sented in  Figure  22.  The  deviation  due  to  thrust  malalignment  is  much 

less  than  that  due  to  mass  unbalance  for  the  values  of  u0  and  OL  used. 

2 y 

The  coupling  of  mass  unbalance  and  thrust  malalignment  is  evident  in  the 
curve  for  Run  5.  The  change  in  yaw  angle  is  much  greater  than  that  in  6 
for  the  parameters  used  as  is  evidenced  by  Figures  22,  23,  and  24. 

A combination  of  variables  often  used  in  deciding  on  the  ‘'goodness1' 
of  a launch  is  the  square  root  of  the  sum  of  the  squares  of  A Q and  i) r. 

This  combination  is  shown  as  a function  of  time  in  Figure  25.  It  is 
merely  fortuitous  that  the  thrust  malalignment  detent  combination 

resulted  in  the  smallest  value  of  ^/(A 0)2  + \j/2  a t the  end  of  the  simula- 
tion time,  since  both  mass  unbalance  and  thrust  malalignment  are  normally 
random  quantities. 

The  variable  6 ^ is  given  as  a function  of  time  in  Figure  26.  Since 
6 ^ is  the  pitch  angle  of  the  rocket  with  respect  to  the  launcher,  the 
time  history  of  6 ^ differs  in  all  cases  (Runs)  from  that  for  A0;  i.e., 

6 ^ includes  the  launcher  pitching  motion  (6^  — Q - C^) • 

Figure  27  illustrates  the  climb  rate  and  lateral  velocity  of  the 

rocket's  center  of  mass.  Since  the  initial  elevation  angle  of  the 

rocket  O'  4-  QL  » 0.093422  radian,  the  rocket's  center  of  mass  will 
20  2e  ’ 

ascend  (this  implies  the  Z becomes  more  negative) , although  gravity  and 
other  factors  will  slow  its  ascent  slightly. 

The  last  figure  (Figure  28)  illustrates  the  fact  that  the  rocket's 
center  of  mass  drops  below  the  initial  aim  axis  due  to  launcher's 
pitching  down,  gravity,  and  the  other  effects  included  in  the  various 
runs  . 
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Figure  13.  Launcher  roll  angle 
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Figure  14.  Launcher  pitch  angle 
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Figure  15.  Launcher  yaw  angle 
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Figure  16.  Body-fixed  angular  velocity  component,  fL  versus  roll  angle, 
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Figure  17.  Body- fixed  angular  velocity  component,  versus  roll  angle 
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Figure  18.  Body- fixed  angular  velocity  components  an<^  ^ 
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Figure  19.  Pitch  rate  versus  time 
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Figure  20.  Yaw  rate  versus  time 
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Figure  21.  Pitch  rate  versus  yaw  rate. 
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Figure  22.  Pitch  angle  versus  yaw  angle. 
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Figure  23.  Rocket  pitch  angle  with  respect  to  aim  axis  versus  time. 
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Figure  24.  Rocket  yaw  augle  versus  time 
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Figure  25.  Square  root  of  the  sum  of  the  squares  of  A G and  i| ; versus  time. 
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Figure  26.  Pitch  angle  of  the  rocket  with  respect  to  the  launcher  versus  time. 
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0.01  0.00  # 0.01  0.02  0.03 

Y (ft/sec) 

Figure  27.  Inertial  components  of  the  velocity  of  the  rocket1 s 
center  of  mass,  Z versus  Y. 
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Figure  28.  Position  of  rocket's  center  of  mass  with  respect  to  the  aim  axis. 


Section  VI.  SUMMARY  AND  CONCLUSIONS 


a.  Summary 

The  problem  of  simulating  the  dynamical  behavior  of  a 
launcher/rocket  system  during  the  entire  launch,  i.e.,  from  a state  of 
rest  with  the  rocket  on  the  launch,  through  the  detent,  guidance,  tipoff, 
and  free- flight  phases,  has  been  considered.  A physical  model  of  the 
system  which  includes  various  factors  that  contribute  to  mallaunch  of 
free  rockets  has  been  described.  Dynamical  and  kinematical  equations 
developed  on  the  basis  of  the  physical  model  have  been  presented  and, 
where  necessary,  converted  from  vector  to  matrix  form.  The  matrix  equa- 
tions have  been  incorporated  in  a digital  computer  code  for  numerical 
solution.  The  computer  code  has  been  used  to  produce  preliminary  results 
and  these  have  been  presented  and  discussed. 

b . Cone lusio ns 


At  the  time  of  the  writing  of  this  report,  the  computer  code 
had  not  been  operational  for  a sufficient  period  of  time  to  provide  the 
data  required  to  reach  any  major  conclusions  regarding  mallaunch.  The 
code  should  be  useful,  however,  in  future  investigations. 

Some  conclusions  which  are  already  essentially  generally  known  are: 

1)  Dynamic  mass  unbalance  is  a primary  contributor  to  mallaunch 
of  spinning  free  rockets. 

2)  Angular  thrust  malalignment  is  also  an  important  mallaunch 
factor . 

3)  Dynamic  mass  unbalance  and  angular  thrust  malalignment  effects 
may  either  tend  to  cancel  or  to  compliment  each  other. 

The  preliminary  results  obtained  indicate  that  the  computer  code 
has  been  properly  implemented  in  its  present  form.  It  should,  however, 
be  modified  to  include  the  effects  of  aerodynamic  forces  and  moments  and 
of  time  varying  rocket  mass  characteristics  during  the  free- flight  phase. 
Such  work  is  planned  under  an  ARD  Grant  to  the  author.* 

For  plans  of  future  work  in  this  area,  the  reader  may  contact  the 
author  or  Dr,  James  L.  Batson,  US  Army  Missile  Command,  ATTN:  AMSMI-RLH, 

Redstone  Arsenal,  Alabama  35805. 


*A  post-LRCP  grant  has  been  awarded  to  the  author.  The  study  funded  by 
this  grant  will  be  conducted  at  Auburn  University,  Auburn,  Alabama  36830. 
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Appendix 

INPUT/OUTPUT  DESCRIPTION 
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CONTROL  CARDS  AND  SAMPLE  INPUT 


OMVLT. 

ACCT  (PN=COCHRAN,  PBC=7294L10000,CC=7800, JN=0001,0P=A2) 
REQUE  S T (TAPE  2 0 , *PF ) 

REWIND (TAPE 20) 

ATTACH (FILE , LHDYN , ID=LHJEC , CY=2) 

REWIND (FILE) 

FILE. 

REWIND (TAPE 20) 

CATALOG (TAPE20 ,LHXXXJC , ID=LHXXX , CY=5) 

EXIT. 

0000000000000000000000 


DEMONSTRATION  RUN  THRUST  MAL  AND  MASS  UNBAL  (ALY=0.0001,  MU2=0.0001) 
FINAL  TIME  DELT  MIN  DELT  IMPRINT  IPLOT 


0.4500 


0.29 
0.00 
0.00 
MU  2 
,00001 


0 


4800.0 

0.0 

1600.0 

0.00 

0.0 

0.0 


0.00 


0.0 

0.0 


0.0001  0.0001  5000  15 

MISSILE  PRINCIPAL  INERTIA  MATRIX 

0.00  0.00 
70.8  0.00 
0.00  70.8 


MU  3 

0.0 


LAUNCHER  INERTIA  MATRIX 

0.0  1600.0 

16000.0  0.0 

0.0  15200.0 

FORE  SHOE  STIFFNESS  MATRIX 

0.00  0.00 

1000000.0  0.0 

0.0  1000000. 

AFT  SHOE  STIFFNESS  MATRIX 

0.00  0.00 

2000000.0  0.0 

0.0  2000000.0 


0 


10.0 

0.1 

0.00 

0.0 

0.0 


LAUNCHER  NATURAL  FREQUENCIES 
2.5 

LAUNCHER  DAMPING  RATIOS 

0.1 

FORE  SHOE  DAMPING  MATRIX 
0.00 
600.0 
0.0 


2.5 


0.1 


0.0 

0.0 

600.0 


0.00 

0.0 

0.0 


AFT  SHOE  DAMPING  MATRIX 
0.00 
800.0 
0.0 


0.00 

0.0 

800.0 
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M 

SM 

RHOF 

RHOA 

ALY 

ALZ 

SIGMAY 

200.0 

7.76 

0.25 

0.25 

0.0001 

0.00 

0.0 

XLA 

YLA 

ZLA 

MU 

SPINT  1 

SPINT  2 

SPINT3 

-3.75 

0.0 

-2.0 

0.1 

0.0 

450.0 

450.0 

T1 

T2 

T3 

T4 

T5 

T6 

T7 

0.05 

0.065 

0.125 

0.1251 

0.126 

0.146 

0.356 

T9 

FI 

F2 

F3 

F4 

0.386 

0.0 

9350.0 

10500.0 

0.0 

DETENT 

SMU 

1800.0 

0.001 

U 

5.75 

0.0 

0.0 

G 

0.0 

0.0 

32.175 

R 

4.0 

0.0 

-2.0 

LF 

-0.5 

0.0 

0.0 

LFA 

-6.583 

0.0 

ALPHAO 

0.0 

0.0 

XLA1 

XLA2 

0.10 

0.0 

0.25  0.25 

I SHOE 
2 

0000000000000000000000 


SIGMAZ 

0.0 

SPINT  4 

450.0 

T8 

0.386 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 


INPUT  FORMATS 


Format 

Card  No. 

Format 

7A10 

23 

3F25.0 

3F10.0 

24 

8F10.0 

3F25.0 

25 

8F10.0 

3F25.0 

26 

8F10.0 

3F25.0 

27 

8F10.0 

2F10.0 

28 

8F10.0 

3F25.0 

29 

3F20.0 

3F25.0 

30 

3F20.0 

3F25.0 

31 

3F20.0 

3F25.0 

32 

3F20.0 

3F25.0 

33 

3F20.0 

3F25.0 

34 

3F20.0 

3F25.0 

35 

2F10.0 

3F25.0 

36 

12 

3F25.0 

3F25.0 

3F25.0 

3F25.0 

3F25.0 

3F25.0 

3F25.0 

3F25.0 
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CONTROL  CARDS  AND  SAMPLE  INPUT 


OMVLT. 

ACCT (PN -COCHRAN , PBC-7294L10000 ,CC=7800 , JN=0001 , 0P=A2) 
REQUEST (TAPE  20 , *PF) 

REWIND (TAPE 20) 

ATTACH ( FILE , LHDYN , ID=LHJEC , CY= 2) 

REWIND (FILE) 

FILE. 

REWIND (TAPE 20) 

CATALOG (TAPE  20 , LHXXX JC , ID=LHXXX , CY=5 ) 

EXIT. 

0000000000000000000000 


DEMONSTRATION  RUN 

THRUST  MAL  AND  MASS 

UNBAL  (ALY=0. 0001 , MU2=0 

0.4500 

0.0001 

500 

15 

0.29 

0.00 

0.00 

0.00 

70.8 

0.00 

0.00 

0.00 

70.8 

0.00001 

4800.0 

0.0 

0.0 

1600.0 

0.0 

16000.0 

0.0 

1600.0 

0.0 

15200.0 

0.00 

0.00 

0.00 

0.0 

1000000.0 

0.0 

0.0 

0.0 

1000000.0 

0.00 

0.0 

0.00 

0.0 

2000000.0 

0.0 

0.0 

0.0 

2000000.0 

10.0 

2.5 

2.5 

0.1 

0.1 

0.1 

0.00 

0.00 

0.00 

0.0 

600.0 

0.0 

0.0 

0.0 

600.0 

0.00 

0.00 

0.00 

0.0 

800.0 

0.0 

0.0 

0.0 

800.0 

200.0 

7.76 

0.25  0.25 

0.0001  0.00  0.0 

-3.75 

0.0 

-2.0  0.1 

0.0  450.0  450 

0.05 

0.065 

0.125  0.1251 

0.126  0.146  0. 

0.386 

0.0 

9350.0  10500.0 

0.0 

1800.0 

5.75 

0.001 

0.0 

0.0 

0.0 

0.0 

32.175 

4.0 

0.0 

-2.0 

-0.5 

0.0 

0.0 

-6.583 

0.0 

0.0 

0.0 

0.10 

0.0 

0.0 

450.0 

0.386 


0.25 


0.25 


0000000000000000000000 
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MENU  DESCRIPTION 


1 - TIME 

2 — DELTA1  = 5..,  X -Displacement  of  point  a in  launcher  basis  (X  -X  ). 

1 Li  Li  Li 

a a 

o 

3 — DELTA2  = &9,  Y -Displacement  of  point  a in  launcher  basis  (Y  - Y ). 

Z L Li  Li 

a a 

o 

4 — DELTA  = Z -Displacement  of  point  a in  launcher  basis  (Z  - Z ). 

3 L L L 

a a 

o 

5 — THETA(l)  = Rocket  roll  angle  relative  to  launcher  axis, 

6 — THETA(2)  = 6 Rocket  pitch  angle  relative  to  launcher  axes, 

7 — THETA(3)  = 6 Rocket  yaw  angle  relative  to  launcher  axes. 

8 — ALPHA(l)  = 0^,  Launcher  roll  angle. 

9 — ALPHA(2)  = 0^,  Launcher  pitch  angle. 

10  — ALPHA(3)  = C ^ , Launcher  yaw  angle. 

11  - DDELTA1  = & 

12  - DDELTA2  = &2 

13  - DDELTA3  = 63 

14  — OMEGAXM  = , x -component  of  Q rocket  angular  velocity. 

1 m — 

15  — OMEGAYM  = Q0>  y^-component  of  ft  rocket  angular  velocity. 

16  — OMEGAZM  = ft^ , z^-component  of  ft  rocket  angular  velocity. 

17  - OMEGAXL  = u)  , x -component  of  w,  launcher  angular  velocity. 

L Li 

18  — OMEGAYL  = w9,  y -component  of  u),  launcher  angle  velocity. 

Z Li 

19  — OMEGAZL  = , z -component  of  u),  launcher  angular  velocity. 

3 Li 

20  — PHI  = <i>  , Rocket  roll  angle  (absolute). 

21  — THETA  = 0,  Rocket  pitch  angle  (absolute). 

22  — PSI  = \(r,  Rocket  yaw  angle  (absolute)  . 
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1 


23  — x = X.  X-coordinate  of  C . 

m 

24  — y = Y,  Y-coordinate  of  C . 

m 

25  — z = Z,  Z-coordinate  of  C . 

m 

26  — u = U,  x -component  of  rocket* s velocity 

m 

27  — v = V,  y -component  of  rocket's  velocity 

m 

28  — w = W,  z^-component  of  rocket's  velocity 

29  - DPDELTA1  = 

30  - DPDELTA2  = &2 

31  - DPDELTA3  = 

32  - DOMEGAXY  = 0 

33  - DOMEGAYM  = 02 

34  - DOMEGAZM  = ^ 

35  - DOMEGAXL  = u). 

36  - DOMEGAYL  = u>2 

37  - DOMEGAZL  = w3 

38  - DPH1  = 0 

39  - DTHETA  = 6 

40  - DPSI  = ijr 

41  - DX  i X 

42  - DY  = Y 

43  - DZ  = Z 

44  - DU  = U 

45  - DV  = V 

46  - DW  = W 
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47  — DELTA  Z = Displacement  of  in  the  Z' -direction. 


48  — DELDZ  = Z + U sin(a^Q  + G^)  component  velocity  of  normal  to 

the  aim  axis. 

49  - DELTHETHA  = 0 - QL  +ao  + A0 

2o  2e 

50  - SRSSA  = M \|r2 
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SYMBOLS 


A 

Direction  cosine  matrix  (launcher  basis  to  rocket  basis). 

a 

Point  on  rocket  centerline  at  the  aft  shoes. 

B 

Direction  cosine  matrix  (inertial  basis  to  launcher 
basis) . 

Oil 

LM 

Oil 

cd 

on 

Damping  coefficient  matrices. 

'l'lVl 

Centroidal  launcher- fixed  coordinate  system. 

C x y z 
m mm  m 

Centroidal  rocket- fixed  coordinate  system. 

C x y z 
m P P P 

Principal,  centroidal,  rocket- fixed  coordinate  system. 

D 

Direction  cosine  matrix  (centroidal  rocket  basis  to 
principal,  centroidal  (rocket  basis). 

F 

Magnitude  of  the  thrust  force. 

F 

—as 

Force  in  the  aft  springs. 

~Fa 

Frictional  force  at  aft  shoes. 

— Ff 

Frictional  force  at  forward  shoes. 

-fs 

Force  in  the  forward  springs. 

F 

— m 

Force  on  the  rocket. 

F 

""D 

Detent  force  on  the  rocket. 

F 

”mF 

Frictional  force  on  the  rocket. 

F 

mg 

Gravitational  force  on  the  rocket. 

F 

Thrust  on  the  rocket  (in  launcher  basis) . 

-T 

Thrust  (rocket  basis). 
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f 


f 

Point  on  the  centerline  of  the  rocket  at  the 
forward  shoes. 

& 

Local  acceleration  of  gravity. 

ys  A A 

I > J jK 

Inertial  unit  vectors. 

II M 

tr* 

Centroidal  inertia  tensor  or  matrix  of  the 
launcher. 

^l/o 

Inertia  tensor  or  matrix  of  the  launcher  about  0. 

6 

Mil 

Centroidal  inertia  tensor  or  matrix  of  the  rocket . 

I 

=pm 

Principal,  centroidal  inertia  matrix  of  the  rocket 

AAA 

VJiA 

Launcher- fixed  unit  vectors. 

Rocket- fixed  unit  vectors. 

I 

Inertia  matrix  of  the  system  about  0 during  the 
spin  up  and  detent  phases. 

Sa-Kf^L 

Stiffness  matrices. 

4 

Vector  from  point  a to  a point  on  the  line  of 
action  of  the  thrust,  Fm. 

i, 

-fa 

Vector  along  the  centerline  of  the  rocket 
directed  from  point  a to  point  f. 

M 

Mass  of  the  launcher. 

m 

Mass  of  the  rocket. 

OXYZ 

Inertial  coordinate  system  with  origin  at  the 
launcher  pivot  point. 

R 

Vector  from  0 to  C . 

Ll 

* XL  ‘l  + yL  JL 

a a 

+ ZL  \ 
a 

Vector  from  C to  a. 

Li 

T / 
— m/a 

Torque  on  the  rocket  about  a. 
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1 


Ti/c 
' m 

Torque  on  the  rocket  about  C . 

m 

T 

~ o 

External  torque  on  the  system  about  0. 

T 

— s 

Spin  torque. 

-f/c 

' m 

Frictional  torque  on  the  rocket  about  C . 

m 

t 

Time. 

u , v , w 

Rocket- fixed  components  of  the  rocket's  velocity. 

AAA 

u = u.i  + €1  + € k Vector  from  a to  C . 

— 1 m yJm  z m m 


V 

Velocity  of  the  rocket. 

X , Y , Z 

Inertial  coordinates  of  C . 

m 

ai,0!2,0!3 

Greek  Characters 

Eulerian  angles  used  in  defining  the  orientation 

of  the  C x y z system. 

L L L L 

a = (a  a2  a3)T 


a ,a 

y z 

Thrust  malalignment  angles. 

Displacements  of  point  a from  its  undeformed 
position. 

e , e 

y z 

Displacements  of  from  the  rocket  centerline. 

Ct  ,£l  ,£l 
\ ly  lz 

Launcher  damping  ratios. 

e 

Rocket  pitch  angle. 

eve2,e3 

Eulerian  angles  used  in  defining  the  orientation 

of  the  C x y z system  relative  to  the  C x y, z, 
mmmm  LLLL 

system. 

n 

Friction  coefficient  (sliding). 

Spin  bearing  friction  coefficient. 
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0 ,0 
y z 

Linear  thrust  malalignment  parameters. 

4> 

Missile  roll  angle. 

Missile  yaw  angle. 

H = i 4-  j 

— 1 m 2 Jm 

+ Q3 

Rocket  angular  velocity. 

s ■ "l  \ + “2  3L 

Launcher  angular  velocity. 

+ “3  \ 

U)  ,00  ,00 

“L  nL  nL 

x y z 

Natural  frequencies  of  the  launcher. 

< >L 

Launcher . 

< >m 

Rocket . 

< >o 

Initial  value  or  undeformed  value. 

< >e 

Change  needed  to  achieve  equilibrium. 

( )_1 

Inverse . 

( )T 

Transpose. 

Pf»  P2 

Distances  from  the  rocket  centerline  at  points 
f and  a,  respectively,  to  the  rocket/launcher 
interface  (rail,  tube  wall,  etc.) 

Other  Notations 

Used  above  a vector  (3x1  matrix),  say  ft,  to 
denote  the  skew- symmetric  matrix, 

/o  -n3  a2  \ 

2 = n3  o -a 

Vn2  S1  0 J 

o 

Used  above  a vector  to  indicate  that  the  time 
derivatives  of  only  the  components  of  the  vector 
are  to  be  taken. 
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A 


E 


Si 


L,m 


Used  above  a vector  or  scalar  to  denote  total  time 
dif ferentation. 

Denotes  a unit  vector. 

Denotes  a vector  or  column  matrix. 

Denotes  a tensor  or  square  matrix. 

Identity  matrix  (3x3). 


Used  with  a vector  with  a circle  above  it  to  denote 
the  system  in  which  the  component  derivative  is 
taken;  e .g . , 


n 

m 


Q-  i 4-  j 4* 

1 m 2 m 3 


k and  — 
m L 


u). 


1 L 


4- 


UK 


A 


4- 


UK 


A 

k m 

Li 
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